home *** CD-ROM | disk | FTP | other *** search
/ Mac-Source 1994 July / Mac-Source_July_1994.iso / C and C++ / Libraries / stringsearch / bmsource / uf.rev.d2.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-05-06  |  3.1 KB  |  156 lines  |  [TEXT/MPS ]
  1. /*
  2.     search routine generated by gen.
  3.     skip=uf, match=rev (using revr), shift=d2
  4. */
  5. /*
  6.  * The authors of this software are Andrew Hume and Daniel Sunday.
  7.  * 
  8.  * Copyright (c) 1991 by AT&T and Daniel Sunday.
  9.  * 
  10.  * Permission to use, copy, modify, and distribute this software for any
  11.  * purpose without fee is hereby granted, provided that this entire notice
  12.  * is included in all copies of any software which is or includes a copy
  13.  * or modification of this software and in all copies of the supporting
  14.  * documentation for such software.
  15.  * 
  16.  * THIS SOFTWARE IS BEING PROVIDED "AS IS", WITHOUT ANY EXPRESS OR IMPLIED
  17.  * WARRANTY.  IN PARTICULAR, NEITHER THE AUTHORS NOR AT&T MAKE ANY
  18.  * REPRESENTATION OR WARRANTY OF ANY KIND CONCERNING THE MERCHANTABILITY
  19.  * OF THIS SOFTWARE OR ITS FITNESS FOR ANY PARTICULAR PURPOSE.
  20.  */
  21.  
  22. #ifndef    CHARTYPE
  23. #define    CHARTYPE    unsigned char
  24. #endif
  25. #define    MAXPAT    256
  26.  
  27. #include    "stats.h"
  28.  
  29. #ifndef    TABTYPE
  30. #define    TABTYPE    long
  31. #endif
  32. typedef TABTYPE Tab;
  33.  
  34. static struct
  35. {
  36.     int patlen;
  37.     CHARTYPE pat[MAXPAT];
  38.     Tab delta[256];
  39.     Tab delta2[257];
  40. } pat;
  41.  
  42. prep(base, m)
  43.     CHARTYPE *base;
  44.     register m;
  45. {
  46.     CHARTYPE *skipc;
  47.     register CHARTYPE *pe, *p;
  48.     register int j;
  49.     register Tab *d;
  50.     register Tab *d2;
  51.     register q1, tp, t, qp, jp, kp;
  52.     Tab f[256], f1[256];
  53.  
  54.     pat.patlen = m;
  55.     if(m > MAXPAT)
  56.         abort();
  57.     memcpy(pat.pat, base, m);
  58.     skipc = 0;
  59.     stats.len = m;
  60.     d = pat.delta;
  61.     for(j = 0; j < 256; j++)
  62.         d[j] = pat.patlen;
  63.     for(p = pat.pat, pe = p+m-1; p < pe; p++)
  64.         d[*p] = pe-p;
  65.     d[*p] = 0;
  66.     skipc = (CHARTYPE *)p;
  67.     d2 = pat.delta2;
  68.     for(j = 1; j < m; j++)
  69.         d2[j] = 2*m-j;
  70.     for(j = m, t = m+1; j > 0; j--, t--){
  71.         f[j] = t;
  72.         while((t <= m) && (base[t-1] != base[j-1])){
  73.             if((m-j) < d2[t])
  74.                 d2[t] = m-j;
  75.             t = f[t];
  76.         }
  77.     }
  78.     q1 = t;
  79.     t = m+1-q1;
  80.     qp = 1;
  81.     for(jp = 1, kp = 0; kp < t; jp++, kp++){
  82.         f1[jp] = kp;
  83.         while((kp >= 1) && (base[jp-1] != base[kp-1]))
  84.             kp = f1[kp];
  85.     }
  86.     while(q1 < m){
  87.         for(j = qp; j <= q1; j++)
  88.             if(m+q1-j < d2[j])
  89.                 d2[j] = m+q1-j;
  90.         qp = q1+1;
  91.         q1 += t-f1[t];
  92.         t = f1[t];
  93.     }
  94. /*for(j=1; j<=m; j++)printf("[%d]=%d ", j, d2[j]); printf("\n");/**/
  95.     d2[0] = m+1;        /* the case where the match succeeded */
  96. }
  97.  
  98. exec(base, n)
  99.     CHARTYPE *base;
  100. {
  101.     int nmatch = 0;
  102.     register CHARTYPE *e, *s;
  103.     register Tab *d0 = pat.delta;
  104.     register k;
  105.     register CHARTYPE *p, *q;
  106.     register CHARTYPE *prev = pat.pat+pat.patlen-1;
  107.     register Tab *d2 = pat.delta2+1;
  108.     register k2;
  109.  
  110.     s = base+pat.patlen-1;
  111.     e = base+n;
  112.     memset(e, pat.pat[pat.patlen-1], pat.patlen);
  113.     while(s < e){
  114. #ifdef    STATS
  115.         k = d0[*s];
  116.         stats.jump++;
  117.         while(k){
  118.             stats.jump++; stats.step[k]++;
  119.             k = d0[*(s += k)];
  120.             stats.jump++; stats.step[k]++;
  121.             k = d0[*(s += k)];
  122.             stats.jump++; stats.step[k]++;
  123.             k = d0[*(s += k)];
  124.         }
  125. #else
  126.         k = d0[*s];
  127.         while(k){
  128.             k = d0[*(s += k)];
  129.             k = d0[*(s += k)];
  130.             k = d0[*(s += k)];
  131.         }
  132. #endif
  133.         if(s >= e)
  134.             break;
  135. #ifdef    STATS
  136.         stats.slow++;
  137. #endif
  138. #define    RH    s
  139.         for(p = prev, q = RH; p > pat.pat; ){
  140. #ifdef    STATS
  141.             stats.cmp++;
  142. #endif
  143.             if(*--q != *--p)
  144.                 goto mismatch;
  145.         }
  146.         nmatch++;
  147.     mismatch:
  148.         k2 = q+d2[p-pat.pat]-RH;
  149. #ifdef    STATS
  150.         stats.step[k2]++; stats.jump++;
  151. #endif
  152.         s += k2;
  153.     }
  154.     return(nmatch);
  155. }
  156.